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ABSTRACT 



Data obtained in the ATIC-2 (Advanced Thin Ionization Calorimeter), 
CREAM (Cosmic Ray Energetics and Mass)) and PAMELA (Payload for Anti- 
matter Matter Exploration and Light-nuclei Astrophysics) experiments suggest 
that elemental interstellar spectra of cosmic rays below the knee at a few times 
10^ GeV are not simple power laws, but they experience hardening at magnetic 
rigidity above about 240 GV. Another essential feature is the difference between 
proton and Helium energy spectra, so that the He/p ratio increases by more than 
50 % in the energy range from 10^ to 10^ GV. We consider the concavity of par- 
ticle spectrum resulting from the nonlinear nature of diffusive shock acceleration 
in supernova remnants (SNR) as a possible reason for the observed spectrum 
hardening. Helium-to-proton ratio increasing with energy can be interpreted as 
a consequence of cosmic ray acceleration by forward and reverse shocks in SNRs. 
The contribution of particles accelerated by reverse shocks makes the concavity 
of the produced overall cosmic ray spectrum more pronounced. The spectra of 
protons and helium nuclei accelerated in SNRs and released into the interstel- 
lar medium are calculated. The derived steady state interstellar spectra are in 
reasonably good agreement with observations. 

Subject headings: acceleration of particles — shock waves — ISM:supernova remnants 
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Introduction 



High-accuracy measurements have revealed deviations of cosmic ray spectra from 
plain po wer laws at energi es 10 to 10 ^ Gev/n. This refers, in particul ar, to the 



ATIC-2 ( 



Panov et al 



( lAdriani et al. 



20091 ) ■ CREAM 



20111 ) experiments; see 



(lAhn et al 



2010 



YoonetaL 



20111 ) and PAMELA 



Lavalld ( 120111 ) for additional references. The general 



conclusion from these measurements is the presence of hardening in the spectra of protons, 
helium and probably heavier nuclei at magnetic rigidity about 240 GV and the increase of 
He/p ratio in the above energy range. 



A number of explanations of these results have been suggested: t 



le hardening may reflec t 



the c ontribution of two distinct populations of cosmic ray sources ( jZatsepin &: Sokolskaya 



20061): the action of numerous sources with the dispersion of source iniection spectra 



([Yuan et al. 



20111 ): fluctuations produced by local cosmic ray sources (iThoudam fc Horandel 



2012); specific conditions of cosmic ray transport in the local bubble (lErlykin fc Wolfendale 



2OI2I ): inte rstellar diffusion with a diffusion coefficient that is not separable in energy 



and space (ITomassetti 



2OI2I ): or the combi nation of diffusion on the cosmic ray induced 



turbulence and the background turbulence (IBlasi et al. 



20121 ). The energy dependent He/p 



ratio may be due to acceleration by shock propagating through the med ium with varyi ng 
chemical composition, including the heliu m wind of a Wolf- Rayet star (iPtuskin et al. 



2010), the stratified material of a bubble ( lOhira fc Iokall201l[ ). or rate of ionization ( iDrury 



201 ll). and the spec ifi cs of thermal ion injec tion in the process of shock acceleration 



flMalkov et al 



2OI2I). 



Biermann et al 



( 120101 ) suggested that the spectral hardening and 
enrichment in heavy nuclei is due to the contribution of a polar cap cosmic-ray component 
produced by supernova explosions in winds of massive progenitor stars. The consideration 
of various interpr e tation s of cosmic ray spectral peculiarities can be found in the paper 



Vladimirov et al. 



(120121 ). 
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In the present work we further dev elop our model of cosmic ray acceleration in 



supernova remnants (jPtuskin et al. 



|2010[ ) (Paper I) to explain both of the required features. 
To a good approximation, the model in Paper I explains the observed overall spectrum of 
cosmic rays in a wide range of energies up to abo ut 10^ GeV. An important im provement 



of the code used in Paper I was made in our work (IZirakashvili &: Ptuskin]l201ll ) (Paper II), 
where particle acceleration by a reverse (backward) shock moving through the material of 
supernova ejecta was included in the calculations. This is in addition to the acceleration by 
the forward shock moving through the circums t ellar r aedium studied in Paper I. The c ode 



described in Paper II was used by 



Zirakashvilil ( 120111 ) ; IZirakashvili fc Aharoniaru (120111 ) for 



modelling particle acceleration in SNRs RXJ1713. 7-3946 and Gas A. 

It is worth noting that similar numerical models of nonlinear production of cosmic rays 
in SNRs, where time- dependent cosmic-ray transport equations are solved together with 
gas-dynamic equatioii s in sperical symrnetry w ere d eveloped by tw'o othe r groups of authors 



( iBerezhko et al. 



1994 : 



Berezhko fc Volk 



20071 ) and (IKang fc Jones 



20061 ). They give results 



analogous to our model for particle acceleration by forward SNR shocks, but they do not 
include acceleration by reverse shocks. Cosmic-ray acceleration by reverse SNR shocks was 



earlier considered by_ 



Telezhinsky et al 



Berezinskii fc PtuskinI ( 1l989l ). The process was recently analysed by 



( I2OI2I ) in a test particle approximation. 



Two effects considered below in the frameworks of our model of particle acceleration 
in SNRs are of particular importance: the nonlinear shock modification that leads to the 
concave spectrum of accelerated particles with a pronounced hardening; and taken into 
account acceleration of supernova ejecta material poor in hydrogen by a reverse shock that 
leads to the difference between the overall proton and helium spectra and increases their 
concavity. 
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2. Modelling of cosmic-ray acceleration in supernova remnants 

Cosmic ray acceleration in shell SNRs proceeds through the diffusive shock acceleration 
mechanism that is a version of the first order Fermi acceleration in the shock vicinity, where 
the gas with a frozen magnetic field is compressing. The fast, charged background particles 
are scattered by random magnetohydrodyna mic (MHD) waves and i nhom ogeneities, and 



they gain energy crossing the shock, see e.g. 



Malkov fc O'C Drurvl (|200l|) for a review. 



The process of efficient acceleration should be modelled simultaneously with the SNR 
hydrodynamics, because the pressure of accelerated particles modifies the gas flow and the 
shock structure that affects the particle spectrum. Also, the current of accelerated particles 
results in the streaming instability that amplifies the background MHD waves. These waves 
in turn determine the value of particle spatial diffusion. 

Hydrodynamical equations for the gas density p(r, t), gas velocity u{r,t), gas pressure 
Pg{r,t), and the equation for the isotropic part of the cosmic-ray proton momentum 
distribution N{r,t,p) in the spherically symmetrical case take the form (see also Paper II): 



dp 1 <9 2 
dt dr ^' 



du du 1 f dPn dPn , 

-^1^ hr^ + ^ ' (2) 



dt dr p \ dr dr J ^ 



dPg dPg ^gPgdr^U dP, 



dN _ 1 d ^2^^^ ^ ^^dN ^dN ^ dN p dr'^w 
^ J.2 Qj, ' ^ dr dr dp Sr^ dr 

r]fS{p-pf) 



Anpim 



-p{Rf + 0, t){Rf - u{Rf + 0, t))6{r - Rf{t)) 
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+ ., 2 PiRb - 0,t)(M(i?, - 0,t) - Rb)5ir - R,it)). (4) 



Here Pc = 4:7t J p'^dpvpN/3 is the cosmic-ray pressure, w{r,t) is the advective velocity 
of cosmic rays, 7^ is the adiabatic index of the gas, and D{r,t,p) is the cosmic-ray diffusion 
coefficient. It was assumed that the diffusive streaming of cosmic rays results in the 
generation of MHD waves. Cosmic-ray particles are scattered by these waves. That is why 
the cosmic-ray advective velocity w may differ from the gas velocity u. In our modelling 
they differ on the value of the radial component of the Alfven velocity calculated in the 
isotropic random magnetic field: w = u + ^aVa/V^, Va = B/y/Anp . Here the factor 
^A describes the possible deviation of the cosmic ray drift velocity from the gas velocity. 
The values = 1 and C,a = ~^ are assumed upstream of the forward and reverse shocks 
respectively, where Alfven waves are generated by the cosmic ray streaming instability 
and propagate in the co rresponding directions. Th e damping of these waves heats the gas 



upstream of the shocks. 



McKenzie fc Voelkl (119821 ). which is described by the last term 



in Eq. (3) and limits the total compression ratios by a number close to 6. In addition, 
we use the value = — 1 downstream of the shocks because the cosmic ray gradient is 
positive in this region. The Alfven drift here strongly influences the slope of the particle 
momentum spectrum, since magnetic fields are compressed downstream of the shocks, and 
the gas velocity in the shock frame falls below the sound speed. This leads, in particular, 
to noticeable steepening of the particle spectrum at the forward shock, see Paper I. 

The two last terms in Eq. (4) correspond to the injection of thermal protons with 
momenta p = Pf, p = Pb and mass m at the fronts of the forward and reverse shocks 
at r = Rf{t) and r = Rb{t), respectively. The indexes / and b are used for quantities 
corresponding to the forward and reverse shock respectively. The dimensionless parameters 
1]^ and 1]^ determine the injection efficiency. The magnetic energy density is relatively small, 
so it does not appear explicitly in Eq. (2). 
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Shocked ejecta and interstellar gas are separated by the contact discontinuity at 
r = Re- The spatial dependence of the magnetic field at r > i?c is taken in the form 

IJn pRf 



B{r,t) = J — ^Jl + {MAVAo/Rfr, r>R,., (5) 
where po is the gas density of the circumstellar mediu m, Vao is the Alfy en velocity there, 



(120051 ) in analysis of 



and Ma is some constant. We employ the results of IVolk et al. 
X-ray radiation from young SNRs, and we assume that magnetic energy density -B^/Svr 
downstream of the shock is 3.5% of the ram pressure poRf that determines the constant 
Ma = 23. According to Eq.(5) the far upstream energy density of the amplified magnetic 



field is a small 0.5/M^ part of the ram pressure. Note that this re^ 



mode lling of the cosmic ray streaming instability in young SNRs ( IZirakashvili fc Ptuskin 



ation is consistent wit 



n 



20081 ). We assume that the magnetic field is compressed in accordance with plasma density 
upstream and downstream of the forward shock. We also assume that there is no decay 
of the magnetic field in the downstream region. The magnetic amplification is weak in 
old remnants with the forward shock speed Rf < MaVao- The magnetic field strength 
B is close to the interstellar value Bq in these SNRs (see Eq. (5)). We assume that the 
magnetic field does not depend on radius downstream of the reverse shock at Rb < r < Rc, 
but depends on density similarly to Eq. (5) upstream of the reverse shock, r < Rf,. 

The diffusion coefficient is of the form 

D = kDb. (6) 

Here Db = vpc/{3ZeB) is the so called Bohm diffusion coefficient for particles of charge Ze 
and velocity v (c is the speed of light). The function /?=(! + (MyiVAo/-R/)^)^ at r > Rc 
and K = 3 at r < i?c approximates the dependence of the diffusion coefficient on shock 
velocity. It corresponds to the case when the MHD turbulence in the shock vicinity is 
amplified by cosmic ray streaming instability, which is balanced by the Kolmogorov-type 
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nonlinearity, see iPtuskin fc Zirakashvilil ( 120051 ) and Paper II for details. We assume that 



the diffusion coefficient is close to its Bohm limit in young SNRs when the shock velocity 
exceeds ~ 700 km/s. This slightly underestimates the diffusion coefficient of particles at 
the very end of the energy spectrum. This is because these particles are scattered by 
small-scale magnetic inhomogenuities generated by the non-resonant streaming instability 
and their diffusion coefficient is several times higher than the Bohm diffusion coefficient 



(jZirakashvili &: Ptuskin 



20081 ). The value of D grows with time because the wave generation 



is becoming less efficient as the shock velocity is decreasing. 

The equation for ions is similar to Eq. (4). For ions with mass number A and mass 
M = Am, it is convenient to use the momentum per nucleon p and the normalization of 
the ion spectra Ni to the nucleon number density. Then, the number density of ions rii is 
Hi = AtcA^^ J p^dpNi. The ion pressure Pi = An J p'^dpvpNi/3 is also taken into account in 
the total cosmic ray pressure Pc- 

The numerical procedure we use to solve Eqs. (l)-(4) was described in detail in Paper 
II. A finite-difference method is employed to Eqs. (f )-(4) upstream and downstream of 
the forward and reverse shocks. A non-uniform numerical grid upstream of the shocks at 
r > Rf and r < allows to resolve small scales of hydrodynamical quantities resulting 
from the pressure gradient of low-energy cosmic rays. The gases compressed at the forward 
and reverse shocks are separated by a contact discontinuity dX r = Rc between the shocks. 
An explicit conservative TVD scheme for hydrodynamical equations (l)-(3) and uniform 
spatial grid are used between the shocks. 
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3. Spectra of cosmic rays produced by supernova remnants 

To demonstrate a possible effect of acceleration by reverse sliocks on cosmic-ray 
composition, we make calculations for two classes of supernovae: Type I where Hydrogen 
is absent in the outer layers of ejecta and Type II where Hydrogen strongly prevails. The 
relative rates of th ese types of supernovae in the Galaxy are 0.46 and 0.54 respectively 



( ISmith et al. 



20111 ). The following characteristics of Type la and Type IIP SNRs as 



representative of these two classes of supernovae are accepted in the calculations. 

Type la SNRs have kinetic energy of explosion E = lO^^erg, number density of the 
surrounding interstellar gas n = 0.1 cm"'^, and mass of ejecta Mgj = 1.4M0. Also important 
for accurate calculations is the index k, which describes the power law density profile 
Ps oc of the outer part of the star that freely expands after supernova explosions; k = 7 
for Type la supernova. Type IIP SNRs have parameters E = lO^^erg, n = 0.1 cm^'^, 
Mej = 8Mq, and k = 12. 

FigH] shows the overall spectra of relativistic nuclei produced in Type la and Type IIP 
SNRs during 10^ years of their evolution by forward (solid lines) and reverse (dash lines) 
shocks. It is assumed for illustration that the interstellar medium and the ejecta of Type II 
SNRs consist of protons; the ejecta of Type I SNRs consist of nuclei with the mass number 
A {A > 1) and the charge Z = A/2. The accepted value of the interstellar magnetic field 
is 5/iG. The pressure of accelerated particles modify the profiles of forward and reverse 
shocks that leads to the concave spectra of cosmic rays. The plotted function is defined 
as Q{p) = Attp'^F{p), where F{p) is the distribution of all accelerated nucleons injected in 
the interstellar medium over the SNR lifetime. The total number of accelerated nucleons 
is / Q{p)dp. The procedure for calculating SNR evolution and simultaneous cosmic-ray 
acceleration was described in-depth in Paper II. One can also find there the plots that show 
the calculated profiles of the gas density, velocity and pressure, and the cosmic ray pressure 
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Fig. 1. — Spectra of nucleons produced during SNR lifetime by forward (solid lines) and 
reverse shocks (dash lines) in Type I (black lines) and Type II (gray lines) SNRs as functions 
of momentum per nucleon p. The case is shown when only protons are accelerated by forward 
shocks in the interstellar medium and by the reverse shocks in Type II SNRs whereas only 
nuclei with mass to charge ratio A/Z = 2 are accelerated by the reverse shocks in Type I 
SNRs. 

in evolving SNRs. 

The spectra of particles accelerated at reverse shocks arc harder than the spectra at 
forward shocks, in spite of the same level of the shock modification for both shocks. This 
is because the shocks propagate in the media with different properties. When the reverse 
shock reaches the fiat part of the ejecta density distribution it propagates in the medium 
density which is decreasing with time. That is why the number of freshly injected ions 
is low in comparison with higher energy particles accelerated earlier, thereby leaving to 
spectral hardening. This effect is absent at the forward shock propagating in the medium 
with constant density. 

Particles accelerated in the numerous SNRs are injected into interstellar space, diffuse 
in galactic magnetic fields, interact with interstellar gas, and finally escape through the 
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cosmic ray halo boundaries into intergalactic space, where cosmi c ray density is negh. 



We employ the pla in diffusion model with a ffat cosmic ray halo fiGinzburg fc Ptuskin 



Strong et al. 



gible. 



1976 



20071 ) for calculations of cosmic ray propagation in the Galaxy. The leaky-box 



approximation to the diffusion model can be used 



of proton and helium intensities, see 



Ptuskin et al 



' or our purpose - the determination 



( 120091 ). The cosmic ray intensity 



obeys the relation I oc u^nQ^X^^ + alma)^^, where Xf. is the escape length (the average 
matter thickness traversed by cosmic rays before they exit from the Galaxy), a is the 
nuclear spallation cross section for a given type of relativistic nuclei moving through the 
interstellar gas, and nia is the mean interstellar atom mass. The escape length is determined 
from the relative abundance of secondary nuclei (primarily from the Boron-to-Carbon 
ratio) in cosmic rays. The approximation formula Xg = 19/3^ g/cm^ at R < 3 GV and 



Ptuskin et al. 



fl2009h (here R is the 



X, = 19/33(i?/3GV)-o-6 g/cm2 at > 3 GV was given in 
particle magnetic rigidity; /3 = v/c). According to the last equation, the resulting spectrum 
is steeper than the source spectrum by 0.6 at high enough energies, but the uncertainty in 
the last value is about 0.1, and statistically accurate measurements of the Boron-to-Carbon 
ratio are not available at energies above ~ 30 GeV/nucleon. In the present calculations we 
assume the dependence X^ oc at -R > 3 GV up to rigidity ~ 10^ GV. It corresponds 

to dependence of the cosmic-ray diffusion coefficient on rigidity in the interstellar medium 
DisM oc R^'^ at relativistic energies. 

The interstellar proton and helium cosmic-ray spectra determined in the framework 
of this propagation model, and use the source functions calculated as described above, 
are presented in Figj2l It is assumed that the reverse shocks accelerate only hydrogen 
nuclei in Type II SNRs and only helium nuclei in 65% of all Type I SNRs; the reminding 
35% are Type Ic SNRs, which do not have significant amounts of hydrogen and helium in 
their ejecta. The forward shocks produced by Type I and Type II supernovae propagating 
through the interstellar gas of standard composition with the relative injection efficiency 
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1 102 10" 10= 

E, GeV/n 



Fig. 2. — Calculated interstellar spectra of protons and helium are shown by thick gray lines. 
Data of the BESS, CAPRICE, AMS-01, ATIC-2 a nd CREAM exp eriments corrected for solar 



Lavalld (1201 ll ). where the corresponding 



modulation effects are taken from the review of 
references and detailed description of observations can be found. The PAMELA experiment 
data are taken from the paper lAdriani et al.l fl2011l) and corrected for solar modulation. 

of hydrogen and helium ions at the forward shocks is chosen to fit the observed He/p 
ratio at the 100 GeV/nucleon reference energy. The shock acceleration of background 
particles with different elemental composition leads to the difference in the overall proton 
and helium spectra seen in Figj2l The forward shock in Type Ib/c SNRs (59% of all Type 
I SNRs) moves through the helium wind of the progenitor Wolf-Rayet star, the last stage 
of a progenitor star evolution for these supernovae, that enriches cosmic rays in helium 
at energies above approximately 10^ GeV/n. The possible contribution of this process is 
shown by dash line in Figl2j 

The calculated spectra in Figj2] exhibit hardening at ~ 240 GV and reproduce the 
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energy dependent He-to-proton ratio in fair agreement with observations. Thus, the 
explanation of these features in terms of acceleration of cosmic rays in SNRs by forward 
and reverse modified shocks seems quite possible. 



4. Discussion and conclusion 



The processes of cosmic-ray acceleration by supernova shocks and subsequent diffusion 
in the interstellar magnetic fields at relativistic energy depend on particle Larmor radius 
so the spectra of different ions might be expected to have similar shapes when express ed 



as fu nctions of p article magnetic rigidity. However , the data of A TIC-2 (iPanov et al 



20091), CREAM flAhn et al. 



201G 



Yoonetal 



201 ll ) and PAMELA flAdriani et al. 



20111) 



experiments showed that the He/p ratio is changing with energy at 10 to 10^ GeV/n. 
Also, the spectra of these ions demonstrate deviations from the plain power laws with 
hardening at about 240 GV. Thus, modern statistically accurate measurements confront 
the traditional power-law paradigm for the source spectrum and escape length of Galactic 
cosmic rays below the knee. A number of possible explanations of these experimental 
results were proposed in the papers cited in the Introduction, but a clear picture does not 
yet emerge. 

In the present paper, the interpretation of both features (spectral hardening and 
changing of the He/p ratio) is provided in the framework of our model of cosmic ray 
acceleration by SNR shocks developed in Papers I and II (the first version of this work 



was presented in the conference paper iPtuskin et al. 



(120111 )). The hardening of cosmic ray 



spectrum results mainly from modification of gas flow in the shock precursor by the cosmic 
ray pressure, which shapes the concave energy spectrum of cosmic rays. This effect is well 



known in the theory o 



fll985h : 



Berezhko et al 



diffusi ve shock acceleration, see e.g. early papers 



Ellison fc Eichler 



(119941 ). It manifests itself here in the presence of Alfven drift of 
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accelerated particles and the strong dependence of cosmic-ray diffusion on shock velocity 
as a pronounced hardening at rigidity ~ 240 GV in the overall spectrum of accelerated 
particles produced by SNRs. The account of particle acceleration by a reverse shock moving 
through the supernova ejecta results in the production of an additional component of cosmic 
rays that has a specific hard energy spectrum and is depleted in hydrogen composition. 
This makes hardening of the overall cosmic ray spectrum more pronounced, and it produces 
the difference between hydrogen and helium spectra. 

The results of our calculations of the interstellar spectra of cosmic-ray hydrogen and 
helium shown in Fig. |5] demonstrate good agreement with observations. It is assumed that 
reverse shocks accelerate only protons in Type II SNRs and only helium nuclei in Type I 
SNRs. Forward shocks accelerate interstellar protons and helium nuclei in both cases. 

Clearly more work is needed and planned to establish the adequacy of the suggested 
scenario. This refers both to the detailed analysis of the dispersion properties of SNRs 
(the SN explosion energy, the mass and complex composition of ejecta, the structure of 
circumstellar medium, etc.) and assumptions about the process of particle acceleration 
(the value of M^, the parametrization of diffusion coefficient, etc.). Also, a single 
power-law dependence of cosmic-ray diffusion in the interstellar magnetic fields on rigidity, 
DisM oc vR^'"^ , used in the present calculations in the entire range of rigidities from 10 to 
10^ GV is most probably a significant oversimplification. The explanation of low anisotropy 
of Galac tic cosmic rays seems incompatible with s uch a strong dependence of diffusion on 



rigidity (jPtuskin et al 



2006 



Blasi &: Amato 



2OI2I). 



In principle, the last difficulty could be relieved by increasing from its Ma = 23 
value taken from Paper I and accepted in the present calculations. This would lead to a 
steeper source spectrum and require weaker dependence of interstellar diffusion on rigidity 
to fit the observations of cosmic ray spectra. Another possibility is a possible suppression 



- 15 - 



of cosmic-ray a nisotropy in the L ocal Bubble with enhanced turbulence where the Solar 



system resides, 



Zirakashvilil (120051 ). 
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